AO-A095  129  BATTELLE  COLUMBUS  LABS  OH  F/«  4/1 

PARTIAL  REFLECTION  EXPERIMENT:  OPERATIONAL  AND  OATA  REDUCTION  T— ETCCU) 
OCT  80  M  A  FLOOD  DAA629-76U-0100 

UNCLASSIFIED  ERADCOM/ASL-CR-80-0100-4  NL 


■joe  nit  COW  Ul  \ 0  9  5 1 2  9 


I 

J 


I 

* 


./  PARTIAL  REFLECTION  EXPERIMENT: 
OPERATIONAL  AND  DATA  REDUCTION  TECHNIQUES  . 


\ 


/  / 

OCTttBt-0580 


By 

I  W-  A-/ Flood 


Electrical  Engineering  Department 
North  Carolina  State  University 
Rayleigh,  North  Carolina 


Under  Contract  DAA-.76-D-0100 


work 


Work  Unit 


_ -r  • 

v75T“‘ftmded  under  DMA  Task 

Reaction  Rates  to  Propagation 

of  O-Region  Ion  Density  and  Conductivity 


Contract  Monitor:  Robert  0.  Olsen 


Approved  for  public  release;  distribution  unlimited 


US  Army  Electronics  Research  and  Development  Command 

ATMOSPHERIC  SCIENCES  LABORATORY 

White  Sands  Missile  Range,  NM  88002 

812  17217 


NOTICES 


\ 


Disclaimers 


The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig¬ 
nated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  be  construed  as  official  Government 
indorsement  or  approval  of  commercial  products  or  services 
referenced  herein. 


Disposition 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not 
return  it  to  the  originator. 


SECURITY  CLASSIFICATION  OF  THIS  page  (Whtn  Palm  Enttrtd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (mid  Sublltlt) 

PARTIAL  REFLECTION  EXPERIMENT: 

OPERATIONAL  AND  DATA  REDUCTION  TECHNIQUES 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

R&D  Final  Report 

S.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORO) 

W.  A.  Flood 

8.  CONTRACT  OR  GRANT  NUMBERfflJ 

DAA-76-D-0100 

-76  dcicc 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Electrical  Engineering  Department7 

North  Carolina  State  University 

Rayleigh,  North  Carolina 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

DA  Task  No.  1L161162B53A/SA2 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Electronics  Research 
and  Development  Command 

AdelDhi.  MD  20783 

12.  REPORT  DATE 

October  1980 

13.  NUMBER  OF  PAGES 

30 

14.  MONITORING  AGENCY  NAME  a  ADDRESSflf  dlfltttnt  from  Controlling  Office; 

Atmospheric  Sciences  Laboratory 

White  Sands  Missile  Range,  NM  88002 

15.  SECURITY  CLASS,  (of  thla  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

1  16.  DISTRIBUTION  STATEMENT  (of  thla  Report)  1 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abetted  on  fared  In  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

Robert  0.  Olsen,  Contract  Monitor 


19.  KEY  WORDS  (Continue  on  roveroo  aide  If  neceaaary  end  Identify  by  block  number) 

D-region 

Electron  densities 
Partial  reflection  sounder 
Middle  atmosphere 

20.  ABSTRACT  pCaKftM  era  rarer—  aide  it  nncaeaery  end.  Identify  by  block  number) 

Operational  and  analytical  procedures  for  reducing  D-region  partial  reflection 
data  obtained  at  White  Sands  Missile  Range  (WSMR)  were  reviewed  and  analyzed 
with  the  expectation  that  this  equipment  would  be  used  to  deduce  D-region 
electron  density  profiles  during  the  total  solar  eclipse  of  26  February  1979 
in  Canada.  The  procedures  devised  as  the  result  of  this  study  are  applicable 


DO  U73  EDITION  OF  I  NOV  SS  IS  OBSOLETE  } 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  <■**•"  Dtlt  Enttrtd) 


r 


i 


S1CUWITY  CLASSIFICATION  OF  THIS  PAOEfWhl  DM*  Mmttnd)  _ _  _ _ 

20.  ABSTRACT  (cont) 

to  data  obtained  during  the  eclipse  at  a  Canadian  site  as  well  as  to  data 
obtained  during  longer  term  operation  at  WSMR. 


2 

SECURITY  CLASSIFICATION  OF  THIS  PAGE(H7)*n  Dmtm  Bnfrmd) 


CONTENTS 


1.1  WSMR  and  the  Partial  Reflection  Experiment.. 

1.2  Data  Screening  and  Reduction......... . 

1.3  Alternate  Means  of  Handling  Data  Reduction.. 

1.4  Deduction  of  Electron  Density  Profiles...... 

1.5  Operation  at  an  Additional  Frequency  at  WSMR 

WSMR  Conclusions  and  Recommendations....... . 

Part  II  -  Eclipse  Measurements................... 

Part  III  -  Summary  and  Conclusions . 

REFERENCES . 


i  ' 

- * _  -  - 


1 . 1  WSMR  and  the  Partial  Reflection  Experiment 

The  WSMR  is  an  ideal  site  for  the  partial  reflection 
experiment.  The  site  is  remarkably  free  from  interference; 
the  lowest  altitude  for  D-region  backscatter  is  set  by  the 
sensitivity  of  the  experimental  equipment  and  by  ground 
clutter  and  not  by  uncontrollable  interfering  signals. 
Measurements  have  shown  that  24-hour  round-the-clock  D-region 
measurements  are  possible  at  this  site. 

The  data  reduction  procedures  in  use  at  WSMR  were 
inadequate  on  several  counts.  There  was  no  screening 
procedure  for  saturated  echoes.  Furthermore,  although  the 
receiver  is  non-linear,  correction  for  receiver  non-linearity 
was  made  only  after  the  A/D  converter  output  had  been 
averaged  rather  than  correcting  each  receiver  output  sample 
and  then  averaging.  Moreover,  at  high  values  of  receiver 
attenuation,  quantizing  noise  became  important. 

Finally,  the  method  of  data  reduction  was  suspect. 
Although  the  Belrose  and  Burke  (1964)  formula  was  used,  the 
collisional  frequency  profile  was  incorrect,  the  refractive 
index  formula  was  not  the  best  available  and  the  data  were 
reduced  by  logarithmically  differencing  measured  values  of 
Ax/A0  at  2-kilometer  altitude  intervals.  The  Physical 
Science  Laboratory  (NMSU)  were  aware  of  the  errors  associated 
with  the  differencing  procedure  and  had  suggested  alternate 
methods  of  data  reduction. 
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1 . 2  Data  Screening  and  Reduction 


The  following  suggestions  were  accepted  by  AS1.  and 
implemented  as  data  screening  procedures  in  the  production 
of  Ax/Aq  profiles  as  a  function  of  altitude. 

1  -  Each  sampled  output  was  corrected  for  receiver  non¬ 

linearity  before  entering  the  averaging  process. 

2  -  At  each  altitude  and  for  each  receiver  attenuation 

setting,  the  number  of  samples  with  an  uncorrected 
amplitude  of  less  than  63  (the  saturated  A-D 
converter  output)  was  made  part  of  the  computer 
print  -  out . 

Before  the  data  were  further  processed,  additional 
screening  procedures  were  instituted.  Assuming  that  the 
amplitudes  of  both  the  0  and  X  wave  returns  are  Rayleigh 
distributed,  one  can  calculate,  as  a  function  of  true  mean 
value  of  the  echo,  the  percent  of  the  samples  which  will 
exceed  a  count  of  63  -  the  maximum  output  of  a  six-bit  A-D 
converter.  The  results  are  presented  in  Table  1  below. 


Table  1 


Percent 

of  samples  exceeding  63 

counts  as 

a 

function  of 

true  mean 

value 

(R.  ) 

v  true' 

ft 

true 

20 

25 

30 

35 

40 

45 

50 

%>  6  3 

.04 

.68 

3.13 

7.85 

14.25 

21.44 

28.75 

A 

*1 

19.972 

24.530 

27.772 

29.212 

29.086 

27.922 

26.290 

*2 

19.998 

24 .960 

29.74  5 

34.157 

38.064 

41.432 

44  .  400 

Also  shown  in  Table  1  are  the  results  of  two  different 
procedures  for  estimating  the  mean  value.  These  last  two 
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entries  will  be  discussed  later.  It  is  clear  that  if  satu¬ 
ration  is  to  affect  less  than  5%  of  the  total 
number  of  measurements,  mean  values  in  excess  of  32 
linearized  counts  are  to  be  eschewed.  The  KSMR  receiver 
used  6  dB  programmed  gain  steps  and  at  altitudes  where  the 
mean  value  was  greater  than  30,  the  average  of  the  next 
attenuation  step  (greater  than  15  counts)  was  used. 

An  examination  of  experimental  data  indicated,  not 
surprisingly,  that  mean  linearized  averages  of  less  than 
10  units  suffered  from  appreciable  quantizing  errors; 
consequently,  these  data  were  not  used.  Good  average  values 
must  be  between  10  and  30  linearized  A/D  values  with  the 
additional  constraint  that  the  number  of  echoes  used  in  the 
average  ought  to  be  greater  than  95%  of  the  total  number  of 
valid  transmitted  pulses. 

The  definition  of  a  valid  transmitted  pulse  was  based 
on  the  fact  that  for  the  normal  D-region  above  KSMR, 
appreciable  echo  amplitudes  were  not  expected  at  altitudes 
as  low  as  50  km.  If  the  echo  amplitude  at  50  km  exceeded 
5  or  10  units,  the  returns  at  all  altitudes  for  that 
transmitted  pulse  were  ruled  invalid  due  to  interference. 
The  low  interference  levels  at  WSMR  were  reflected  in  the 
fact  that  Ax/A0  averages  computed  from  valid  pulse  levels 
of  5  and  then  10  units  were  virtually  identical. 

The  final  data  screening  procedures  applied  to  the 
linearized  A/D  outputs  were 

1)  The  mean  value  had  to  lie  between  10  and  30  counts. 


2)  The  mean  value  had  to  be  computed  from  at  least 
95 o  of  valid  transmitted  pulses. 

3)  Individual  amplitudes  greater  than  62  counts 
were  discarded. 

These  screening  criteria  were  complementary  for  the 
most  part.  Average  amplitudes  above  30  counts  almost 
invariably  had  more  than  51  of  the  echoes  rejected  for 
saturation.  Average  echo  amplitudes  less  than  10  counts 
differed  significantly  from  the  averages  in  the  next  6  dB 
gain  column  showing  the  effects  of  quantization  noise. 

The  screening  procedures  were  readily  imp  1 ement ab 1 e  and 
produced  consistent  AXM0  profiles  from  the  raw  data.  The 
A^/Ao  profiles  so  produced  differed  from  similar  profiles 
obtained  by  the  author  at  Raleigh,  North  Carolina,  in  that 
the  WSMR  curves  maximized  at  78  kilometers  whereas  the 
Raleigh  curves  maximized  between  70  and  72  kilometers.  It 
was  suggested  that  receiver  delay  was  not  properly  accounted 
for  in  the  WSMR  data.  Tests,  conducted  by  PSL  personnel, 
confirmed  that  this  was  the  case  and  when  receiver  delay 
corrections  were  incorporated  into  the  data  reduction 
procedures,  A^/Aq  profiles,  similar  to  those  obtained  at 
Raleigh,  resulted. 

1 . 3  Alternate  Means  of  Handling  Data  Reduction 

In  the  data  reduction  process  outlined  above,  all 
echoes  whose  instantaneous  amplitude  exceeded  62  were  dis¬ 
carded  and  not  used  in  the  averaging  process.  This  should 
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result  in  a  bias  in  the  average  to  lower  than  true 
values.  We  propose  to  evaluate  the  magnitude  of  this  bias. 
The  third  row  in  Table  1  indicates  the  estimate  of  the 
average  as  a  function  of  the  true  average  value  ^true- 
It  can  be  seen  that  so  long  as  less  than  5%  of  the  echoes 
exceed  62  (part  of  our  screening  procedure)  the  errors  in 
estimating  the  true  mean  value  are  less  than  7 %  — -  an 
acceptable  value.  However,  Dr.  Dave  Mott  of  PSL  suggested 
that  we  keep  a  count  of  the  number  of  echoes  which  exceed 
62  counts  and,  instead  of  discarding  them,  count  them  as  63. 
The  estimate  of  the  true  mean  value  for  this  case  is  denoted 
by  R2  and  is  shown  in  the  fourth  row  of  Table  1.  It  is 
clear  that  so  long  as  less  than  15%  of  the  echoes  are 
saturated,  this  technique  results  in  less  bias  in  the 
estimate  of  the  average  over  a  wider  range  of  true  mean 
values.  I  recommend  that  this  procedure  be  adopted  in 
future  data  reduction  procedures  although  I  would  recommend 
that  apparent  mean  values  greater  than  35  be  used  with 
caution  —  the  percentage  of  saturated  echoes  will  be 
greater  than  8  -  15%  and  possible  combinations  of  sampling 
errors  and  saturated  echoes  may  result  in  spurious  results. 

1 . 4  Deduction  of  Electron  Density  Profiles 

The  data  screening  procedures  outlined  in  section  1.2 
result  in  curves  of  Ax/A0  which  must  be  further  processed 
to  produce  electron  density  profiles.  The  experimental 
technique  is  based  upon  a  theoretical  prediction  of  the  ratio 
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of  amplitude  of  the  extraordinary  return  to  the  amplitude 
of  the  ordinary  wave  return.  Theory  predicts  that  the 
measured  ratio  as  a  function  of  altitude  ought  to  behave  as 


|RxCh)  i 2 

■pewr 


exp 


dz } 


(1) 


is  the  mean  squared  amplitude  of  the  X  wave  return 
is  the  mean  squared  amplitude  of  the  0  wave  return 

is  the  ratio  of  the  backs  cat t e r  i ng  cross  section 
(per  unit  volume)  for  X  and  0  waves,  respectively 
and  is  independent  of  electron  density. 

the  imaginary  part  of  the  refractive  index  for 
the  extraordinary  wave. 

the  imaginary  part  of  the  refractive  index  for 
the  ordinary  wave. 

are  the  absorption  indices  for  the  ordinary  and 
extraordinary  polarizations,  respectively. 

|E  (h) |2 

The  ratio  -  is  proportional  to  the  square  of  the 

TE7KTT2 

ratio  of  the  mean  amplitude  of  the  X  wave  (A^)  to  the  mean 

amplitude  of  the  0  wave  (An)  at  that  altitude  and  is  a 

|Rx(h) | * 

measured  quantity.  -  is  a  theoretical  ratio.  The 

TRo~00  1 2 

difference  between  the  measured  and  theoretical  ratio  is  due 
to  the  differential  absorption  of  the  X  and  0  wave  polar¬ 
izations  and  is  a  function  of  the  electron  density  as  a 
function  of  altitude. 

The  formula  in  equation  (1)  is  essentially  that  used 
by  Gardner  and  Pawsey  (1954)  in  their  original  paper 


TEJ2 
I  Rx  I  2 

htt2 
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explicating  the  partial  reflection  technique.  It  is  also  the 
fundamental  formula  proposed  by  Bel  rose  and  Burke  in  their 
1964  paper.  Belrose  and  Burke  were  the  first  to  use  the 
Sen-Wyller  formulation  for  the  refractive  index  of  an 
ionized,  magnetized  plasma  -Gardner  and  Pawsey  had  used  the 
older  Appleton-Hartree  formulation  which  assumed  a  constant 
mean  free  path  independent  of  electron  energy.  Flood  (1968) 
modified  this  basic  formula  to  account  for  absorption  within 
the  scattering  volume.  In  addition  to  this  correction  term, 
Flood  proposed  to  invert  equation  (1)  using  a  least  squared 
fit  to  a  polynomial  solution  for  N(h),  the  electron  density 
profile  as  a  function  of  altitude,  in  place  of  solving  a 
series  of  coupled  difference  equations  for  equation  (1) 
evaluated  at  successive  altitudes.  The  advantage  of  the 
least  squared  polynomial  fit  is  that  it  minimizes  the 
propagation  of  numerical  and  experimental  errors  in  the 
inversion  process  associated  with  deducing  electron  density 
profiles  from  experimental  data. 

The  major  changes  in  the  WSMR  data  reduction  procedures 
instituted  during  this  project  were  primarily  the  use  of  a 
least  squared  polynomial  solution  in  place  of  the  coupled 
difference  equations  previously  used.  Additional  changes 
included  a  revised  electron-neutral  collisional  frequency 
profile  and  a  newer  and  more  accurate  approximation  to  the 
complete  Sen-h’yller  theory.  These  changes  will  be  discussed 
in  order. 
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A  simplified  collisional  frequency  profile  suggested  by 
Davis  was  used  for  data  reduction  (1969).  At  Raleigh  we  had  been 
using  a  3-piece  fit  to  the  1962  COSPAR  atmosphere  and  there  were 
obvious  differences  between  the  two  profiles.  The  Atmospheric 
Sciences  Laboratory  made  available  the  results  of  9  months  of 
Robin  Sphere  data  obtained  at  WSMR.  [The  Robin  Sphere  is  a 
rocket  launched  sphere  with  an  accurately  known  ballistic 
coefficient.  Radar  tracking  of  the  sphere,  released  at  rocket 
apogee,  provides  estimates  of  atmospheric  density,  pressure,  and 
ultimately  electron-neutral  collisional  frequency  as  a  function 
of  altitude.]  The  Robin  Sphere  data  were  analyzed  and  were  found 
to  agree  very  closely  with  the  three-piece  fit  to  the  COSPAR 
atmosphere.  This  collisional  frequency  profile  was  therefore 
incorporated  into  the  WSMR  electron  density  reduction  routine. 

The  Sen-Wyller  formula  for  the  refractive  index  of  an 
ionized  magnetized  gas  does  not  lend  itself  very  well  to  computer 
inversion  of  equation  (1)  and  solutions  for  electron  density  as  a 
function  of  altitude.  Beirose  and  Burke  proposed  a 
q ua s i -1 ongi t ud i na 1  (Q-L)  approximation  which  was  appropriate  for 
the  high  geomagnetic  latitude  associated  with  theit  Ottowa 
location.  Unfortunately,  there  was  no  way  to  logically  extend 
their  formation  to  lower  latitudes.  Flood  (1974  unpublished) 
suggested  the  following  Q-L  approximation  for  D-region  mid-  as 
well  as  h igh- la t i t ude  refractive  indices. 


12 


n  =  1 


sS  |  T2<h>  *  i  522  T4<h> 


,  NChJe'- 

nx  2me  ojv  (h)  |_ 

0  u 


Tjth)  +  j  5/2  T3(h) 


(2a) 

(2b) 


where  N(h)  is  the  electron  density  at  the  altitude  h 


m  =  mass  of  an  electron  =  9.11x10 
e  =  charge  of  electron  =  1.602x10 


31 


19 


kg 

coulombs 


cQ  =  permittivity  of  free  space  = 

8.85x10  ^2  farads/meter 

v  (h)  =  electron  neutral  collisional  frequency 

w  =  angular  operating  frequency  (radians/sec) 


Tl(h)  =  Ayx  C3/2(yx)  ♦  Byo  C3/2(yo) 

+  D  C3y2( y) 

(3a) 

T2(h)  =  Ayo  C3/2<yQ)  +  Byx  C3/2(yx) 

+  D  C3y 2 (y) 

(3b) 

T3(h)  -  A  C5/2(yx)  ♦  B  C5/2(yo)  +  D 

C5/2(y) 

(3c) 

T4(h)  -  A  C5/2(y0)  +  B  C5/2(yx)  ♦  D 

C5/2(y) 

(3d) 

A  =  cos?  ((f>/2)  -  ~  s in 2  ( d>/ 2 ) 

(4a) 

B  =  sin2(4>/2)  -  ^  sin2(4>/2) 

(4b) 

D  =  j  sin2(4>) 

(4c) 

<p  =  acute  angle  between  the  direction  of  propagation 
the  earth's  magnetic  field 

and 

C  4  d) 

U)- (!>„  U  +  10m 

-  H  .  v  -  H  Y  = 

U) 

(4e) 

yx  \T(hy  ’  yo  v  (h)  ’  y 

v(h) 

eB 


=  angular  electron  gyrof requency  =  m 

2 

B  =  magnetic  flux  density  of  earth's  field  (webers/m  :Teslas) 


f  co  P  X 

11  x*  e 

C  (v)  =  semi-conductor  integral  =  — y  -  -  *■  dx 

P  ’  P'{  * 

0 

The  polynomial  approximations  to  the  semiconductor 
integral  (O'Hara  1964)  can  be  used  profitably  in  the  data 

reduction  process. 


(4f) 
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It  now  remains  to  evaluate  the  ratio  of  backseat  t  rn  ng 
cross-sections  for  the  extraordinary  to  ordinary  polarization1' 


In  terms  of  equation  (1) 


I  Rx (h)  I  2  2 

R  (h) 


Tj (h)  ♦  ^  T|(h) 

T|(h)  +  ^  T*(h) 


(5) 


2 

Flood's  1968  correction  would  multiply  R  (h)  by: 

no  (h)  sinh(-y-  n*  ct) 
n*(h)  sinh(-^  n^  ct) 

where  t  is  the  resolvable  system  pulse  length. 


If  the  pulse  length  t  is  short  enough  (less  than  25 
microseconds)  and  the  electron  densities  not  too  high, 

Flood's  correction  term  does  not  differ  significantly  from 
unity.  In  instances  of  unusually  high  absorption,  this 
correction  term  may  be  important  even  for  pulse  lengths  as 
short  as  20-25  microseconds.  Note  that  the  effective  pulse 
length  of  the  system  is  determined  by  the  entire  receiver 
response  and  not  just  by  the  transmitted  pulse  length.  The 
effective  pulse  length  is  given  by  the  convolution  of  the 
impulse  response  of  the  receiver  (a  measurable  quantity)  and 
the  transmitted  pulse  shape  (also  measurable).  The  overall 
response  of  the  receiver  can  be  measured  by  noting  the  3  dB 
width  of  the  returned  pulse  from  a  classical  E-region  return 
wherein  echoes  are  found  only  over  a  narrow  range  of  alti¬ 
tudes.  At  any  rate,  it  appeared  that  for  the  normal  daytime 
D-region  above  WSMR,  the  Flood  correction  is  not  necessary 
over  the  altitude  region  60-90  kilometers,  using  20  microsecom 
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effective  pulse  lengths.  The  current  procedure  in  use  at 
WSMR  neglects  Flood's  correction  term.  The  form  of  equation 
(1)  in  use  at  WSMR  is  therefore  given  by 


T 1  (h)  +~3T 


T  ,  (h) 


T22(h)  + 


2  S 

T 


Tj(h) 


r  8tt 

exp  {-  — 


(nJ-iV)  dzl  ...(7) 


where  all  symbols  have  been  defined  in  equations  2,3, 
4  ,  S  and  6 . 


is  known  from  measurements  at  2-kilometer  altitude 
increments.  T^  ,  T2,  T^,  T^  are  known  functions  of  altitude, 
electron  gyrof requency ,  operational  frequency  and  magnetic 
dip  angle.  The  differential  absorption  term, 

n  fh 

exp  {-  -r-  I  (n1-!!1)  dz),  is  the  only  term  involving  electron 
A  j  X  O 


density. 

There  are  two  ways  of  inverting  (7).  The  first,  used 

by  Gardner  and  Pawsey  (1953)  and  Belrose  and  Burke  (1964),  is 

to  take  the  logarithm  of  equation  (7)  and  then  obtain  a 

series  of  coupled  difference  equations  which  can  be  solved 

for  a  slab-wise  approximation  to  the  electron  density  at 

each  height  h.  Flood,  in  1968,  suggested  use  of  a  polynomial 

solution  for  elect ron  dens i ty  to  provide  a  least  squared  fit 

Ax.2 

to  the  logarithm  of  (- — >  .  The  least  squared  fitting 

Ao 

procedure  can  be  applied  to  equation  (7)  whether  Flood's 
correction  term  is  used  or  not.  Furthermore,  one  has  the 
option  of  weighting  the  contributions  to  the  mean  squared 
residuals.  In  particular,  Flood  (1968)  suggested  minimizing 
the  mean  squared  percentage  error.  The  current  WSMR  procedure 
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uses  a  polynomial  approximation  for  the  electron  density 

profile  with  uniform  weighting  of  the  mean  square  errors. 

The  advantage  of  the  polynomial  solution  is  that  each  measured 

point  is  recognized  for  what  it  is  —  an  estimate  of  the 
Ax  2 

true  value  of  ( — )  at  that  altitude  rather  than  forcing  a 

Ao 

solution  to  pass  through  every  experimental  point  —  some  of 
which  may  he  significantly  in  error.  The  coupled  difference 
equation  solution  suggested  by  Gardner  and  Pausey  is  very 
prone  to  numerical  instability.  A  slight  error  at  a  low 
altitude  leads  to  large  errors  (negative  electron  densities) 
at  higher  altitudes.  The  polynomial  least  squared  fit,  when 
the  order  of  the  polynomial  is  restricted  to  less  than  half 

Ax 

the  number  of  heights  at  which  —  is  determined,  does  not 

A 

suffer  from  this  problem. 


1 . 5  Operation  at  an  Additional  Frequency  at  WSMR 

The  equipment  installed  at  WSMR  consists  of  two  separate 
transmitters  capable  of  200  kw  peak  pulse  power  outputs,  two 
receivers  and  sufficient  digital  data  recording  apparatus  to 
enable  partial  reflection  operation  at  two  different  frequen¬ 
cies.  The  advantage  of  such  an  arrangement  is  that  a  lower 
frequency  (say,  2.66  MHz)  can  be  used  to  deduce  electron 
densities  over  the  range,  say,  of  60  -  88  kilometers  and  a 
higher  frequency  can  be  used  to  explore  electron  densities 
over  the  range  (say)  of  75  -  95  kilometers.  Candidate 
frequencies  for  the  second  system  were  4.5  and  6.0  MHz. 
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Calculations  of  the  expected  —  ratio  expected  at  4.5 

A0 

and  6.0  MHz  operation  at  WSMR  using  average  summer  and 
winter  D-region  electron  density  profiles  (as  experienced 
in  Raleigh,  NC)  were  made  and  the  results  presented  in 
Figures  1  and  2.  It  is  apparent  that  if  meaningful  data  are 
to  be  collected  on  a  diurnal  basis,  the  second  operating 
frequency  should  be  between  2.66  and  4.5  MHz. 


ALTITUDE  IN  KILOMETERS 


Figure  2.  Calculated  ratios  of  A  /A  and  R  /R  . 

X  O  X  o 
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WSMR  Conclusions  and  Recommendat i ons 


The  data  screening  procedures  described  above  can 
minimize  the  task  of  deducing  valid  electron  density  piolile- 
from  the  raw  data.  There  is  one  additional  constraint  to 
the  overall  data  reduction  problem  which  needs  to  he  addressed 
the  order  of  the  polynomial  in  electron  density  used  to 

Ax  2 

describe  the )  data  points.  It  is  clear  that  in  the 
o 

D-region,  the  electron  density  increases  very  roughly, 

exponentially  with  altitude.  We  require  at  least  seven  valid 

A  2 

contiguous  measurements  of(^p-)  as  a  function  of  altitude  for 

o 

the  data  run  to  be  considered  for  further  analysis.  If  there 

are  at  least  seven  valid  contiguous  altitude  measurements  of 
A 

the  order  of  the  polynomial  for  electron  density  is 
o 

constrained  to  be  less  than  (n-l)/2  where  n  is  the  number 
A 

of  contiguous  measurements.  Although  the  constraints 
o 

enumerated  above  may  reject  a  large  quantity  of  data 
collected  at  a  normal  D-region  station,  WSMR,  because  of  its 
isolation,  can  in  fact  collect  a  large  quantity  of  valid  data. 

If  WSMR  institutes  a  second  operational  frequency,  then 
a  frequency  of  approximately  3.5  MHz  is  strongly  recommended. 
At  WSMR,  these  two  frequencies  should  be  able  to  provide 
reasonable  noon-time  electron  density  profiles  over  the  region 
66  to  94  kilometers.  Furthermore,  substantially  more  valid 

data  could  be  procured  if  an  eight-bit  A/D  converter  were  sub¬ 
stituted  for  the  six-bit  converter  currently  in  use.  It  is 
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recommended  that  future  programs  at  WSMR  utilize  a  gain- 
programmable,  "super- 1 inear"  receiver  of  the  type  in  use 
at  Raleigh.  The  major  portions  of  a  gain-programmable, 
"super- 1 inear”  receiver  have  been  delivered  to  PSL. 


Part  II  -  Eclipse  Measurements 

2.1  The  Physical  Sciences  Laboratory,  New  Mexico  State 
University,  was  to  operate  the  partial  reflection  experiment 
in  Canada  during  the  total  solar  eclipse  of  February  26,  1979. 
In  anticipation  of  problems  which  might  be  encountered  in 
reducing  partial  reflection  data  from  the  eclipse  site, 
equipment  calibration  and  data  reduction  procedures  in  use 
at  the  WSMR  field  site  were  reviewed  and  the  results  of  this 
review  have  been  presented  as  Part  I  of  this  report. 

It  was  not  expected  that  the  P .  R-.  experiment  would 
provide  ground-based  observations  of  the  rapid  changes  in 
electron  density  at  altitudes  between  65  and  80  kilometers 
as  a  result  of  the  obscuration  of  the  sun  during  the  eclipse. 
The  P.R.  experiment  was  expected  to  provide  continuous  (as 
contrasted  to  the  instantaneous  rocket  measurements) 
measurements  of  electron  density  above  80  kilometers  and 
continuous,  "filtered"  measurements  of  electron  density 
below  80  kilometers.  Furthermore,  although  a  few  rockets 
can  be  fired  on  "control"  days  to  ascertain  normal  variations 
of  electron  densities  in  the  D-region,  the  P.R.  experiment 
offered  the  opportunity  to  continuously  measure  D-region 
electron  densities  at  almost  all  times  after  the  equipment 
has  been  installed. 

The  D-region  results  from  the  Eclipse  Experiment,  while 
not  unexpected  because  of  the  high  geomagnetic  latitude  of 
the  experiment  site,  were  significantly  different  from  those 
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results  expected  from  an  eclipse  experiment  at  a  mid-  or 
low-latitude  site.  These  differences  cause  some  difficulty 
in  data  analysis.  An  example  of  the  kind  of  difficulty 
can  be  seen  from  an  inspection  of  Table  2. 

Table  2  lists  the  average  value  of  the  A-D  output 

(after  linearization)  for  ordinary  and  extraordinary  returns 

as  a  function  of  apparent  altitude.  This  is  a  1-minute 

sample  taken  at  1044,  February  26,  1979  —  approximately 

9  minutes  before  second  contact.  One  hundred  pulses  for  each 

polarization  and  each  attenuation  were  used  to  compute  the 

averages.  An  asterisk  in  the  corner  of  each  listing  indicates 

that  more  than  5%  of  the  echo  amplitudes  were  saturated  (count  = 

6^  units).  These  values  were  not  used  in  the  calculation  of  the 

~  ratios. 

A 

o 

If  we  restrict  our  analysis  to  those  altitudes  at 
which  less  than  51  of  the  echoes  were  saturated  and  in  which 
the  mean  values  were  between  10  and  30  units,  then  we  can 
see  that  the  6  dB  incremental  receiver  gain  change  is 
reflected  in  the  data.  For  example,  using  ordinary  wave 
polarization,  the  first  altitude  at  which  the  eqho  amplitude 
at  both  6  dB  and  12  dB  attenuation  levels  are  both  within 
the  range  10  -  30  linearized  counts  is  a  66  km.  The  average 
amplitude  at  the  6  dB  attenuation  level  is  23.128  counts. 

The  corresponding  average  amplitude  for  the  12  dB  attenuation 
level  is  11.571  counts  —  differing  by  only  .061  from  the 
expected  value  of  11.564.  Furthermore,  look  at  the  AQ 
values  in  the  12  and  18  dB  attenuation  level  columns. 
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TABLE  2.  CORRECTED  MEAN  RF.C FIVER  OUTPUT  AS  A  FUNCTION  OF  ALTITUDE  POLARIZATION  AND 
ATTENUATION  SETTING 
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A.  ffiihifftiin 


ES 


If  the  count  in  the  18  dB  column  is  less  than  10,  it  is  very 
difficult  to  see  a  factor  of  2  (6  dB  difference)  in  the  two 
column  values;  note,  however,  that  at  74  and  80  kilometer? 
where  both  counts  are  within  the  10-30  count  range,  the 
significance  of  the  6  dB  (factor  of  2)  gain  change  is  clearly 
shown.  The  6  dB  and  12  dB  columns  illustrate  even  more 
clearly  that  if  both  values  lie  between  10  and  30  counts, 
it  is  only  sampling  error  which  disturbs  the  factor  of 
2  expected  difference.  We  take  these  observations  as  further 
evidence  of  the  desirability  of  confining  the  valid  entries 
to  those  values  which  lie  between  10  and  30  counts  (all  data 

having  less  than  51  saturated  values) . 

A 

The  ratioswere  computed  only  for  those  entries 
o 

which  met  the  screening  criteria  above.  When  values  were 
available  for  two  attenuation  levels  —  such  as  56  km 
for  the  ordinary  wave  in  Table  2  —  the  two  values  were 
averaged  (giving  a  6  dB  value  of  23.135  for  example  in 
Table  2 ,  66  km) . 

In  looking  over  the  eclipse  data  furnished  by  ASL,  the 
author  noted  that  at  no  altitudes  (using  say  6  dB  ordinary 
wave  values  in  Table  2)  did  the  tabulated  values  go 
significantly  below  a  count  of  10  units.  One  might  be 
tempted  to  say  that  for  this  gain  setting  something  like  a 
count  of  10  constitutes  a  "noise"  level  such  that  all  6  dB 
values  listed  ought  to  have  a  count  of  10  subtracted  in 
order  to  account  for  the  "noise"  contribution.  [For  the 
12  dB  column,  the  least  value  listed  is  approximately  3.5 
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and  3.5  might  be  considered  the  "noise"  contribution  to  the 

A 

12  dB  values.]  One  can  construct  new  profiles  corrc.  ted 

Ao 

for  noise  in  this  fashion  and  this  author  lias  done  so 

even  though  at  Raleigh  we  never  used  a  "noise"  correction. 

Upon  reflection,  at  least  for  the  eclipse  data,  such  a 

correction  is  not  warranted.  In  fact,  after  making  such 

"noise"  corrections,  the  factor  of  2  between  succersive  6, 

12,  18  dB  columns  disappears  at  essentially  all  altitudes. 

The  preliminary  electron  density  profiles  using  a  noise 

correction  produced  by  this  author  in  meetings  with  ASL 

personnel,  cannot  be  supported  and  should  be  ignored. 

It  is  clear  that  1-minute  averages  of  echo  amplitudes 

suffer  from  sampling  errors  but  that  longer  averaging  times 

will  prevent  the  recognition  of  any  rapid  ec 1 i pse -produced 

changes  in  electron  density.  Consequently,  it  is  recommended 

A 

that  1-minute  averages  of  produced  by  the  screening 

o 

procedures  already  discussed  (less  than  5%  saturation  and 

counts  between  10  and  30)  should  in  turn  be  averaged  with  a 

running  weighted  average:  a  3-minute  average  with  a  1,2,1 

weighting  readily  suggests  itself.  These  weighted  -r— 

Ao 

averages  should  then  be  subjected  to  a  polynomial  inversion 
and  in  view  of  the  high  absorption  experienced  at  the 
eclipse  site,  the  inclusion  of  the  "Flood  sinh  term"  is 
recommended.  The  polynomial  solution  can  then  be  weighted 

Ax 

to  minimize  the  percentage  error  in  the  j—  profile. 
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A  . 

The  profile  of  —■  available  from  Table  2  is  interesting 
o 

(and  in  some  ways  inexplicable)  from  several  points  of  view. 

A  . 

At  a  normal  mid-latitude  station  one  expects  the  profile 

o 

to  increase  monoton i cal ly  to  a  maximum  value  and  thereafter 

decrease.  Below  56  kilometers  and  above  78  kilometers  this 

expected  behavior  is  not  reflected  in  the  experimental  data. 

Measurements  of  the  cone  angle  of  arrival  of  D-region  echoes 

at  Raleigh,  NC,  show  that  when  this  expected  pattern  is  not 

followed,  oblique  echoes  and  altitude  smearing  are  frequently 

encountered.  Provision  to  measure  the  cone  angle  of  arrival 

of  the  D-region  echoes  was  not  available  at  the  eclipse  site 

A 

so  that  the  cause  of  the  inexplicable  behavior  of  the 

o 

ratio  above  78  kilometers  cannot  be  documented  at  this  time. 


I  would  not  include  experimental  values  outside  the  58-to 

78-kilometer  region  were  I  reducing  the  data  shown  in  Table  2. 

Unless  there  were  extraordinarily  strong  effects  for  the 

differential  absorption  within  the  scattering  volume  — 

Flood's  correction  —  it  is  difficult  to  explain  the  non- 

A 

monotonic  behavior  of  the  profile  above  78  kilometers 

o 

and  blind  data  reduction  which  attempted  to  produce  electron 
density  profiles  from  these  data  would  be  thoroughly  suspect. 

In  summary,  the  data  from  the  eclipse  of  2/26/79, 
taken  from  a  site  in  Canada,  show  strong  particle  ionization 
characteristics  —  so  strong  that  the  solar  eclipse  may  have 
had  only  a  minor  effect  upon  the  D-region  electron  densities 
during  the  course  of  the  eclipse.  The  "control  day"  data  for 


the  most  part  also  show  strong  particle  ionization  effects. 
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While  the  P.R.  data  gathered  during  the  eclipse  expe¬ 
dition  may  not  have  fulfilled  the  expectations  of  a  solar 
controlled  experiment,  the  P.R.  data,  backed  up  by  a  most 
unique  array  of  in-situ  rocket  measurements,  can  very  well 
provide  the  best  documented  and  measured  example  of  a  parti¬ 
cle  precipitation  event  yet  extant.  The  eclipse  experiments 
include  rocket  borne  measurements  of  all  kinds  of  input  flux 
to  the  D-region,  P.R.  electron  density  profiles,  rocket  borne 
measurements  of  electron  density  and  rocket  borne  mass 
spectrometry  measurements  of  positive  and  negative  ion  density. 

Tn  many  ways  the  eclipse  data  may  be  more  definitive 
in  checking  the  DAIRCHEM  code's  predictions  than  the 
morphology  of  a  classical  solar  eclipse  might  have  been. 

The  data  collected  during  the  eclipse  expedition  are 
appropriate  to  a  particle  ionized  D-region  —  something 
closer  to  the  real  requirement  for  a  DAIRCHEM  code  check  than 
could  be  provided  by  a  simple  solar  eclipse. 
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Part  III 


-  Summary  and  Conclusions 

A  quick  look  at  the  eclipse  P.R.  data  implies  that  the 
eclipsed  D-region  was  dominated  by  exceptionally  strong 
particle  ionization  sources.  Even  pre-eclipse  and  post- 

A 

eclipse  runs  frequently  show  non-solar  controlled  profiles. 

^o 

The  P.R.  results  coupled  with  the  rocket-borne  in-situ 
measurements  of  the  eclipse  experiment  in  Canada  can  very 
well  provide  more  interesting  scientific  data  and  better 
tests  of  the  DAIRCHEM  code  predictions  than  could  a  con¬ 
ventional  lower  latitude  eclipse  experiment. 

It  is  recommended,  in  view  of  the  high  absorption,  that 
the  Flood  "sinh"  correction  term  be  incorporated  into  the 
analysis.  It  is  also  recommended  that  no  "noise"  correction 
be  made  to  the  data.  Finally,  following  the  suggestion  of 
Dave  Mott  to  include  saturated  data  points  (at  a  63  value) 
provided  less  than,  say,  15%  of  the  samples  at  that  altitude 
are  saturated  is  endorsed.  The  interference  situation  at 
the  experimental  site  in  Canada  was  such  that  the  imposition 
of  valid  data  screening  levels  at  50  kilometers  seems  not 
necessary,  and  misleading  in  that  there  were  real  D-region 
echoes  at  50  kilometers. 
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